Spin-density- wave instabilities in the organic conductor (TMTSF)2C104: Role of anion 

ordering. 
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We study the spin-density-wave (SDW) instabilities in the quasi-one-dimensional conductor 
(TMTSF)2C104. The orientational order of the anions CIO4 doubles the unit cell and leads to 
the presence of two electronic bands at the Fermi level. From the Ginzburg-Landau expansion of 
the free energy, we determine the low-temperature phase diagram as a function of the strength of 
the Coulomb potential due to the anions. Upon increasing the anion potential, we first find a SDW 
phase corresponding to an interband pairing. This SDW phase is rapidly suppressed, the metallic 
phase being then stable down to zero temperature. The SDW instability is restored when the anion 
potential becomes of the order of the interchain hopping amplitude. The metal-SDW transition 
corresponds to an intraband pairing which leaves half of the Fermi surface metallic. At lower tem- 
perature, a second transition, corresponding to the other intraband pairing, takes place and opens 
a gap on the whole Fermi surface. We discuss the consequence of our results for the experimental 
phase diagram of (TMTSF)2C104 at high magnetic field. 

PACS Numbers: 75.30Fv, 74.70. Kn,81.30.Dz 



I. INTRODUCTION 



The organic conductors of the Bcchgaard salt family 
(TMTSF)2X (where TMTSF stands for tetramethylte- 
traselenafulvalene and X=PF6, CIO4...) exhibit a very 
rich phase diagram when temperature, magnetic field or 
pressure are varied.EJ One of the most remarkable phe- 
nomena is the existence of a series of spin-density-wave 
(SDW) phaseSjin presence of a moderate magnetic field 
of a few TcslaJj These phases are separated by first-order 
transitions and exhibit a quantization of the Hall effect: 
aj,y — ~2Ne^/h per layer of TMTSF molecules, where 
the integer N varies at each phase transition. 

Accorjding to the so-called quantizcd-nesting model 
(QNM),Erlj the formation of the magnetic- field-induced 
spin-density-wave (FISDW) phases results from an in- 
terplay between the nesting properties of the quasi-one- 
dimensional (QID) Fermi surface and the quantization 
of the electronic orbits in magnetic field. Although the 
QNM explains the quantization of the Hall effcctQEl and 
most features of the phase diagram, deviations from the 
theoretical predictions have been observed at liiigh mag- 
netic field in the compound (TMTSF)2C104Ja"i In the 
last FISDW phase, when the magnetic field exceeds 18 
Tcsla, the second-order metal-SDW transition which oc- 
curs at Tc ~ 5.5 K is followed by a SDW-SDW tran- 
sition at ~ 3.5 K. It is believed that the existence 
of this low-temperature SDW phase is due to the ori- 
entational ordering of the (non-centrosymmetric) anions 
CIO4 which occurs at Tao — 24 K in slowly cooled (re- 
laxed) samples.Elia Nevertheless, to ouij-knowlcdge, there 
is no satisfying theoretical descriptiora of the phase di- 
agram oij-(TfMTSF)2C104 at high field in spite of recent 
progress. E3t13 



In this paper, we study the effect of anion ordering 
on the SDW phase in the absence of a magnetic field. 
In relaxed samples, the ground-state of (TMTSF)2C104 
is superconducting in the absence of a magnetic field. 
Therefore, our study of the SDW in presence of anion 
ordering does not apply directly to the experimental sit- 
uation. Nevertheless, it is the first step towards the un- 
derstanding of the behavior of (TMTSF)2C10^|at[±|iffih| 
field. Our work shows that in previous analysialOliil't^^t^ 
(with or without magnetic field) , the anion ordering was 
not accounted for properly. Because of the anion po- 
tential, the SDW order parameter necessarily has two 
Fourier components. As a result, the metal-SDW tran- 
sition temperature is determined by a generalized Stoner 
criterion. This point has been systematically overlooked, 
and this calls for a revision of previous works. 

The results obtained in this paper are based on a sim- 
ple model, where the anion ordering is assumed to create 
an electrostatic potential V {—V) on even (odd) chains. 
Moreover, the crystal structure is taken to be orthorhom- 
bic, while the actual structure of the Bechgaard salts is 
triclinic. In the spirit of the QNM, we expect such a sim- 
ple model (with a few unknown parameters) to correctly 
describe the physics of (TMTSF)2C104. The anion po- 
tential ztV doubles the crystal periodicity in the trans- 
verse direction. This leads to a reduced Brillouin zone 
with two electronic bands crossing the Fermi level (Fig. |l|) 
in qualitative agreement with the actual Fermi surface 
of (TMTSF[)aC104 as obtained from quantum chemistry 
calculation.li3 Instead of a single (best) nesting vector, 
there are three possible nesting vectors: Qinter (inter- 
band pairing), Q-|_ and Q_ (intraband pairing). The 
instability which does occur at low temperature depends 
on the ratio V/tb between the anion potential and the 
interchain hopping amplitude tb- It has recently been 
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shown, both-from quantum chemistry calculationlla and 
experimental^, that V can be of the order of ti,. 

In the next section, we calculate the electron-hole 
susceptibility within the random-phase approximation 
(RPA) and thus obtain the transition temperature be- 
tween the metallic phase and the SDW phase. The 
complete phase diagram is obtained from the Ginzburg- 
Landau expansion of the free energy (Sec. HI). For a 
weak anion potential {V <C h), we find that the SDW 
phase corresponds to the interband pairing (Qinter)- The 
transition temperature is strongly suppressed by the an- 
ion potential and in general vanishes above a critical 
value of V. The SDW instability is restored when V be- 
comes of the order of the interchain hopping amplitude 
{V ^ tb). The metal-SDW transition corresponds to the 
intraband pairing Q+ (or Q_). Half of the Fermi surface 
remains gapless, so that the SDW phase is metallic. At 
lower temperature, a second SDW instability occurs at 
the (intraband) nesting vector Q_ (or Q+), thus opening 
a gap on the whole Fermi surface. 

Some of oufi .conclusions agree with the results of 
Kishigi et a/.t2rt3. In particular, these authors have 
shown that the anion potential may stabilize SDW phases 
with wave vectors Qinter, Q+ (or Q_ and Q+).Ea 



II. INSTABILITY OF THE METALLIC PHASE 

In QID materials, the Fermi surface consists of two 
slightly warped open sheets. As a result, in the vicinity 
of the Fermi level, the electron dispersion is well approx- 
imated as 



f^aik^, ky) = vpiaka: - kp) + t^{kyh), 



(2.1) 



where k^ and ky are the electron momenta along and 
across the conducting chains, and h is the interchain 
spacing. Here and in the rest of the paper, we ne- 
glect the third direction {z axis) which does not play 
an important role for our purpose, and use natural units 
Ti = kB = c = 1. In Eq. (pH]), the longitudinal elec- 
tron dispersion is linearized in k^ in the vicinity of the 
two ID Fermi points ifc^^, and vp — 2ataSm(kpa) is 
the corresponding Fermi velocity. 4 is the transfer inte- 
gral along the chain and a the lattice spacing, a = — , + 
corresponds to the left and the right Fermi sheets. The 
transverse dispersion t±{kyb) is given by 

t^ikyb) = tf^ikyh) + tT'\kyb) (2.2) 

tT'^ikyb) ^ -2*26 cos{2kyb) (2.3) 
t°^'^{kyb) = -24 cos{kyb) - 2^^ cos{3kyb) (2.4) 

where tnb is the transfer integral for electron hopping 
to the nth neighboring chain. t°^'^ and t°^°'^ corre- 
spond to odd and even n, respectively. For a simple 
model with only nearest-neighbor hopping, i.e. when 
t2b = tsb = 0, the linearized dispersion (2.1) satisfies 
the property e-(k) = — e+(k + Qo), which corresponds 



to a perfect nesting of the Fermi surface at wave vec- 
tor Qo = (2fci?,7r/5). For t2b, Hb 7^ 0, the nesting be- 
comes imperfect. Consequently, the nesting vector shifts 
to Q = [2kF + Sqx, n/b + 6qy) and the SDW phase occurs 
with lower transition temperature. Generally, t^b is ne- 
glected since it is very small, i.e., t^b/t2b ^ 1- However, 
as we shall see, plays an important role in the SDW 
phases with intraband nesting vectors (5+ and Q-, and 
should therefore be retained. 

The anion potential in QID systems can be most sim- 
ply modeled as 14nion = V'(— 1)", where n is the chain 
index. In the presence of the anion potential, the Hamil- 
tonian of the system in the absence of electron-electron 
interaction can be written as 



i?o= / d^r^'ij)i^{v)\vF{akx-kF) + t^{kyb) 

+V^(-l)"lV'a.(r), (2.5) 



where the V'qo-'s are fermionic operators for right [a = -I-) 
and left (a — — ) moving particles, and a =t, i is the spin 
index, r = {x,nb), and kx,ky are momentum operators 
along X and y. 



The Hamiltonian (2.5) can be diagonalized to obtain 
the eigenf unctions and energy eigenvalues: 



ipK-r/2 



vpiak^ - kp) + 4 ' 



(2.6) 
(2.7) 



where j — ±, A ~ L^Ly is the area of the system, 
K = (0,7r/fe), and ky e] — 7r/2&, 7r/26] refiecting period 
doubling along y due to the presence of the anion poten- 
tial. The transverse momentum dependent part of the 
energy is given by 

< ^ J^]V^ + [tf''{kyb)Y +tT'\kyb). (2.8) 

The presence of the anion potential splits the energy dis- 
persions into two bands with energies and e~ as 

shown in Fig. |l|. The factors ^■'^ are given by 



1/2 



it 



7^- = 1 1 



1/2 



1/2 



tf^ikyb) 



(2.9) 



The amplitude of the wavefunction at r = (x, nb) de- 
pends on the factors ^jj' and is given by 



iek«i = i7r + (-ir7ri. 



(2.10) 
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Electronic states in the + (— ) band have a higher prob- ^ 1 j ^ |__|_ 

abihty amphtude on even (odd) chains. This locahza- Xaa- (li 1 + I^) — ~aZ^ \ \Xaa (l^i 1, 
tion becomes very important when V ^ v^if, (assum- 
ing ^ if,, which is obviously the case for a realistic x {^t'^it + 7zt^ it^ 

. \ y y y ^y "■y 9 

dispersion law). In this regime, even and odd chains 
tend to decouple, and the dispersion law reduces to 

Using Eqs. (2.6) and (p??]), we obtain the Green func- 
tions in the absence of electron-electron interaction: 




where 



Gaa(r,r;tj„) = }_^— , (2.11) , r^ST ( ■ o Y (■ / 



where a;„ is a fermionic Matsubara frequency. (2-16) 



A. Bare susceptibilities 

Due to the presence of two electronic bands, there are 
three possible spin-density-wave (SDW) instabilities with 
wave- vectors Q+, Q_. and Qintor as shown in Fig. |^. 
Q+ and Q_ are intraband nesting vectors which satisfy 
{Qy)± ~ 7i'/2&. Qintcr is the inter-band nesting vector 
(with ((3y)intor ~ ■"'Z^)- In the absence of the anion po- 
tential, the SDW instability occurs with Q = Qintor- To 
obtain a quantitative description of these possible insta- 
bilities and the phase diagram for the system, we calcu- 
late the susceptibility x(q, q') in the particle-hole channel 
within the RPA. First, let us consider the susceptibilities 
in the absence of electron-electron interaction: 

r'; r - r') = (T.A„.(r, T)AL(r', t'))o, (2.12) 

where Kaa{Y,T) = ?/'t_(r, r)-0Q„(r, t), t is an imaginary 
time, and is th e time-ordering operator. The mean 
value in Eq. (2.12) has to be taken with the Hamilto- 
nian Ho. We use the notation a = —a, and a =i,t for 
a =t,i- In frequency doma in, y *^ can be expressed in 
terms of the Green function ( 2.11 ) as 



xL(i"7i"';Pn) -T^GQ^(r,r';w„)Gaff(r',r;t 



(2.13) 



where p„ is a bosonic Matsubara frequency. For studying 
the instabilities of the metallic phase, it is sufficient to 
compute the static susceptibilities x°„(r,r';p „ = 0) 



xL(r,r'). Using Eqs. (|2j), (^J), (0), and we 
find 

xL(q, q) = ^ E I (1^' + (k, q)] 



+- 



X Wfcy ^ky-qy ^ky~qy^ky 



(2.14) 



Notice that due to the presence of the anion potential, the 
static susceptibilities have a non-zero off-diag onal com- 
ponen t x^{^.^ q + K). The sum over in Eqs. (2.14) and 
15 ) can be analytically calculated using 



bL 



^^X^:(k,q) 



7V(0) 



log 



27^0 



-Re* 



, (2.17) 



but the ky sum needs to be evaluated numerically. In 



Eq. (2.17), N{0) = l/irvph is the density of states per 
spin at Fermi energy, Eq ~ ta an ultraviolet cutoff en- 
ergy, 7 ~ 1.783 the exponential of the Euler constant, 
* the digamma function, and Re ^' means real part of 
Note that when V ^ 0, XQa(k, k + q) = and the 
susceptibility becomes diagonal in momentum space. 

To find the effect of the anion potential V on the bare 
susceptibilities, we plot x°(q, q) and x'^(qi q + K) for dif- 
ferent values of V/tb in Figs. |^^. Here = X+| ~ X+ji 
and q is chosen such that ^ 2k p. For V/tt = 0, x° 
is diagonal in momentum space and the peak of x''(qi q) 
is located at Qinter with (Qinter)y ~ T^/b (Fig. |). The 
position of the peak moves away from ir/b when devi- 
ations from perfect nesting due to t2b and tsb become 
important. As we increase V/tb, x"(q, q + K) becomes 
non-zero and develops peaks at qy w -in: /2b as shown, 
for V/tb = 1, in Fig. ^. The maximum of X^lq^q) ('-S- 
X^(Qintor, Qintcr)) rcduccs in height and two additional 
peaks develop at Q_|- and Q_ (with {Q±)y = ■n/2b). The 
development of these additional peaks in x° (q, q) can 
be seen by comparing Figs. and ^. When V is strong 
enough, the maximum of x moves from Qintcr to Q-|_ 
or Q_. Wc therefore expect the SDW wave- vector to 
shift from Qintcr to Q+ and/or Q_ as V exceeds a criti- 
cal value (assuming that the anion potential V does not 
suppress the SDW instability). Our results for the diag- 
onal susceptibility x"(qiq) are similar to those of Ref. 

m. 
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B. RPA calculation 

To find the critical value of the anion potential and to 
obtain the phase diagram, we now compute the suscepti- 
bilities for the interacting system within RPA. We model 
the interaction-using the g-ology model keeping only (72 
to be non-zero:c3 



f ^ / d^r^LWV^L'W^a.'WV'a.W. (2.18) 



Using the interaction Hamiltonian given by Eq. ( 2.18 ), 
we obtain the susceptibilities by summing the RPA dia- 
grams shown in Fig. |[ 

xa<T(q,q') = xL(q'q') + .92X!^"T(q>q")xaa(q'^q')• 
q" 

(2.19) 



Substituting Eqs. (|2.14D and ( |2.15| ) in Eq. (|2.19| ), we find 
that XacT (q7q') is non-zero only when q' = q or q + K 
and is given by 

xaa(q,q) = {xL(q,q)[i-.92xL(q + K,q + K)] 
+ff2[xL(q,q + K)]2}/z?, (2.20) 

Xa.(q,q + K) = xL(q,q + K)/A (2.21) 

^ = [1 - .92xL(q, q)][i - 52xL(q + k, q + k)] 

-52[xL(q,q + K)]^ (2.22) 
A SDW instability occurs when the susceptibihtics 



,(q,q) and y°^(q, q +K ) diver 



when the de- 
nominator D in Eqs. (|2.20| ) and ( ^.21]) vanishes. This 



yields a generalized Stoner criterion for a SDW instabil- 
ity at wave vector Q: 

[1 - ff2X°(Q, Q)] [1 - ,92X°(Q + K, Q -f K)] 
-52[X°(Q,Q + K)]2 = 0, (2.23) 

where Q is chosen such that Qx ^ 2kp. This equa- 
tion is the same for Q and Q -I- K (reflecting the fact 
that the SDW contains Fourier components at Q and 
Q -|- K). We choose to label the SDW phase by the 
wave vector Q such that x''(Q, Q) > x°(Q + K, Q + K). 
Note that the off-diagonal component of the susceptibil- 
ity, i.e. x''(Q,Q + K), has been systematically over- 
looked in previous works on the effect of anion ordering 
in (TMTSF)2C104. As a result, the simple (but wrong) 
criterion 1 — (72x''(Qi Q) = has been used to determine 
the SDW transition temperature instead of the general- 
ized Stoner criterion given in Eq. (2.2S). 

Eq. (2.23) has to be numerically solved to obtain the 
transition temperature Tc{V) as a function of the an- 
ion potential V. The qualitative nature of the phase di- 
agram, however, can be understood qualitatively from 
Eq. ( p.23| ) without numerical computation. For V = 0, 
the off-diagonal susceptibility x°(Qj Q + K) = and Eq. 



(|^ reduces to 1 - g2X°(Q,Q) = 0. The SDW insta- 
bility occurs with the nesting wave-vector Qintcr corre- 
sponding to the maximum of the bare susceptibility x*^- 
When V increases, the wave functions become more and 
more localized on even or odd chains. Since Qintcr cor- 
responds to pairing of electron and hole on chains of 
opposite parity, Tc drops, until at some critical value 
of the anion potential, Vd, there is no SDW instabil- 
ity with Q = Qintcr- [Note that the anion potential 
docs not affect the nesting at Qintcr, which is limited 
by t2b in our model.] Simultaneously, the pairing with 
Q = Q+ and/or Q_ becomes increasingly favorable, due 
to both the wave function localization and an improved 
nesting, and eventually at a value of the anion poten- 
tial, Vc2 ^ tb, the SDW phase with Q = Q+ and/or Q„ 
wins over the metallic phase. Indeed, when V 3> tb, 
the even and the odd chains tend to decouple. The 
dispersion law of the two electronic bands (-1- and — ) 
becomes ^ — '<^F{ctkx — kp) + jV + t'^^'^'^{kyb), and 

T,{V) r]°^ - (27^o/7r)exp(-2/7V(0)g2), where ri°) 
is the transition temperature for a system with perfect 
nesting. Higher-order harmonics in the transverse dis- 
persion law, such as t4b (not considered in this paper), 
would introduce deviations from perfect nesting that sur- 
vive when V 00. They are however expected to be very 
small and can be safely discarded. If Vc2 > Ki, there will 
be a metallic region in the phase diagram with no SDW 
instability for Vd < V < Vc2', otherwise, there will be a 
transition from the SDW phase with Q = Qintcr to the 
SDW phase with Q = Q+ and/or Q_. 

Our numerical calculations show that both the scenar- 
ios are possible. The phase diagram with t2b = ^36 = is 
shown in Fig. |^. In this case, there is a transition between 
SDW phases with Q = Qintcr and Q = Q± without any 
intermediate metallic phase. Moreover, the SDW's with 
Q = Q-i- have the same transition temperature. The 
transition between the SDW phases with Qintcr and Q± 
occurs a.t V — Vc ^ l.ltb, as indicated by the thick ver- 
tical line in Fig. ^. 

The phase diagram with more realistic parameters, 
shown in Fig. is quite different. With t2b = O.ltf, 
and t^b = 0, the SDW phase with Qintcr persists up to 
Vci ~ O.Atb. For V > Vci, we find a metallic phase till 
V ^ Vc2 ^ tb- For V > Vc2.. the SDW phase with Q± 
is stabilized. When t2b is strong enough, the SDW at 
Qintcr is Suppressed {i.e. Vd = 0), but the transition 
temperature of the intraband SDW's is not affected. In 
this case, the metallic phase is stable at T = up to Vc2 
where the SDW with wave vector Q+ or Q_ sets in. 

The degeneracy between Q-|_ and Q_ is lifted by the 
presence of a non-zero t^b (Fig. 0). For t^b > 0, the SDW 
with wave vector Q_ has a higher transition tempera- 
ture, i.e. T~ > . The opposite is true when < 0. 
It is clear that when V ^ tb, the two SDW's will coexist 
at low temperature. In this regime, even and odd chains 
are essentially decoupled, and the two (intraband) in- 
stabilities take place almost independently at transition 
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temperatures T+ and . This leads to a phase with 
two coexisting SDW's below T™^'' ~ T+ < T~ (assum- 
ing % > 0). 

Fig. I shows the ratio J = (1 + tan(6'))/(l - tan(6')) 
of t he SD W amplitudes on even and odd chains (see 
Sec. |[IIDD . 



In the presence of two SDW's, the order parameter 
AQ,cr(r) = (AQCT(r)) takes the general form 



A„„(r) 



(3.1) 



Note that the SDW with wave vector Qj corresponds 
to a spin modulation with Fourier components Qj and 
+ K. The order parameter ( p.l[ ), which corresponds 
to electron-hole pairs with opposite spins, assumes the 
SDW's to be polarized in the {x,y) plane, i.e. {Sz{r)) = 
In the preceding section, we have shown that two SDW 0. The spin mo dulati on (S(r)) is discussed in more de- 



III. GINZBURG-LANDAU EXPANSION 



instabilities (at wave vectors Q+ and Q_) may occur tail below (Sec. |III D| ). The relation Aaff(r) = A*^(r) 
when V is strong enough. In order to study in more implies that the complex order parameters A^^ satisfy 
detail the possibility of coexisting SDW's, we derive the Ag*^ = A^p. . 



Ginzburg-Landau expansion of the free energy and then 
deduce the phase diagram. 



J 



Up to quartic order in the order parameter, the free 
energy (per unit surface) is given by 



^ = 92^1 I ^I^^WI'+ffl^E / rf'r-idV2A;(ri)A„(r2)G„T(ri,r2,c.)G5Kr2,ri,c.) 
+■92^ E / rf'^irf'^2dV3d2r4A:(ri)A„(r2)A:(r3)A„(r4) 



xG'Q|(ri,r2, w)Gai(r2,r3,w)G'at(r3,r4, w)Gax(r4,ri,a;), 
where A£j(r) = Ac-f(r). 



(3.2) 



A. Quadratic contribution F2 

Let us first consider the quadratic contribution F-i to the free energy. One easily obtains 

-ff2X°(Q,,Q, +K)\ / Ai+ 
, +K,Q, +K); Ar 



W - n \^(A3+* AJ-*\f l-32x"(Qj,Qj) -927 
F,-g,^{A^ ,A„ )(,-g2X°(Q,+K,Q,) l-.g2X°(Q, 



Introducing the new order parameters Uja,Vja defined by 



Ar 



cos{9j) —sm{9j)\ f Uja 



sm 



cos 



with 9. J e] - 7r/4, 7r/4] and 



we obtain the diagonal form 



tan(26' 



2x"(Q„Q,+K) 



X°(Q.+K,Q,+K)-x°(Q„Q,)' 

^2 = E(^^i"^"i'+vi^j"i') 



where 



'± - ^^(2 - .g2X°(Q„ Q,) - 92X\Q, + K, Q, + K)) ± | 

: ([X°(Q,, Q,) - X°(Q, + K, Q, + K)]2 + 4[x°(Q„ Q, + K)] 

r 



A± = f (2 - .92X°(Q„ Q,) - 92X\Q, + K, Q, + K)) ± ^sgn(x°(Q, + K, Q, + K) - x°(Q„ Q,)) 

2\l/2 



(3.3) 



(3.4) 



(3.5) 



(3.6) 



(3.7) 
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The transition temperature is determined by 

minAj*(T^) = 0. (3.8) 

^From Eqs. ( |3.7| ) and ( |3.8| ), one easily recovers the gen- 
eralized Stoner criterion obtained in the preceding sec- 
tion [Eq. (2.23)]. Since we hav e ass umed X°(QijQj) > 
X°(Qj + K,Qj + K) (see Sec. |IB|), minA^^ = A+ and 
the order parameter of the transition is Uj^ ■ For T < T^!, 
we then have Vja = 0, i.e. 



Ar =tan(0,)Ai+. 



(3.9) 



This equation determines the relative amplitude of th e 
spin modulation on even and odd chains (see Sec. HI D| ). 
If > T+ {i.e. t3b > 0), then u-a -f- when T <T-. 
To determine whether we can have also u+a 7^ at lower 
temperature (i.e. coexistence of two SDW's), we must 
analyze the quartic contribution to the free energy. 



of a SDW at wave vector Q^, the nesting is good (and 
becomes better with increasing V). Diagrams contribut- 
ing to C mix the bands -I- and — . A typical diagram is 
shown in Fig. Contrary to diagrams of Fig. it is 
not possible to have all electronic states near the Fermi 
surface and at least one of them has an energy of order 
V with respect to the Fermi level. For instance, in the 
diagram of Fig. ^d, the state (k— Qj + Q_j , j') has an en- 
ergy of order V with respect to the Fermi level when the 
states (k, j) and (k— Qj, j) he near the Fermi surface. As 
a result, these diagrams turn out to be of order {T/V)'^ 
with respect to those contributing to Bj (see Appendix 
^). As shown in the next section, this ensures that at 
low temperature the two order parameters u+q and u_q, 
coexist. 



C. Phase diagram 



B. Quartic contribution F4, 

The calculation of the quartic contribution is some- 
what lengthy, and we only give the main results. Mor e 
details can be f ound in Appendix ^ Using Eqs. (2.6), 
(2.11) and ( |3.lD , we rewrite the quartic part of the free 
energy [Eq. (p.S|)] as 



Fa 



E 



J1---J4,P1---P4 

'XBa{jl,j2,j3,j4;Pl,P2,P3,P4) 

y,/\^jlPl* /\^j2P2 /\^j3P3* /\^j4.P4 



(3.10) 



The coefficients Ba are defined in Appendix We now 
express as a function of the order parameters Uja and 
Vja. Since the formation of the SDW at wave vector Qj 
corresponds to the order parameter Uja , we take Vja = 
in the following. We then obtain a free energy of the form 



+ C|u+„U_a|^|. 



(3.11) 



The expression of the coefficients Bj and C is given in 
Appendix The i mpor tant point here is that the inter- 
action term in Eq. ( 3.1l| ) is weak, i.e. 




(3.12) 



This result is easily obtained by considering the diagrams 
that contribute to Bj or C. The main contribution to Bj 
comes from the type of diagrams shown in Fig. Qj be- 
ing the best nesting vector for the pairing in the band j, 
the value of this diagram is essentially determined by the 
nesting properties. Since the quadratic term F2 predicts 
an instability of the metallic phase against the formation 



Collecting the results from the preceding sections, we 
obtain the following free energy: 

a j 

(3.13) 

where Aj = A~^. 

The free energy F is analyzed in Appendix ^ Assum- 
ing that T~ > T+ {i.e. t^h > 0), we find a second-order 
metal-SDW transition at temperature T~ ^ below which 
U-a 7^ and u+a = 0. This transition opens a gap on the 
— band, while the + band remains gapless. This SDW 
phase is metallic. Note that the minimum of F corre- 
sponds to = \uj-\. This im plies that the SDW is 
linearly polarized (see Sec. Ill D ). 

If AB^B^ — > 0, a second transition takes place at 



Tt + i{Tt-T; 



(3.14) 



where 7 is a constant of order T"^ /V"^ . Below T^°°^, both 
u-a and u+Q are finite, m+q opens a gap on the + band, 
making the whole Fermi surface gapped, so that this low- 
temperature SDW is truly insulating. The transition is 
also of second order, since the order parameters u_q. and 
M+Q vary continuously at the transition. 

Since C/Bj = 0{T^ /V^), the condition 4B+B_ - 
> is satisfied. Therefore, there is always coexis- 
tence of two SDW' s at l ow temperature when the pairing 
is intraband. Eq. (3.14) shows that T^°^^ is well approx- 
imated by T+. 

SDW phases with wave vector Q_, Q+, or Q_ 
and- Q J. 1 1 have been previously obtained by Kishigi et 
a/.E£rt5ll3 However, the overall phase diagram (and in 
particular the existence of a SDW-metal-SDW transi- 
tion) as a function of V has not been derived before. 



6 



D. Spin modulation (S(r)) 

In this section, wc determine the spin modulation 
(S(r)) in the different SDW phases of the phase dia- 
gram. S(r) = (l/2)X;a,a.^' V'L(i")T^a'V^a<T'(i-) is the 
spin-density operator and t — {Tx,Ty,Tz) stands for the 

PauH matrices. Noting that S- = Sx — iSy = Aqi, 
wc obtain 

a,j,p 

= Ai+e^'^^'- [1 + (-1)" tan(0,)] , (3.15) 

Oi,j 

where we have used Vja = to obtain the second hnc of 
Eq. ( ^.15| ). We consider the more general case where two 
SDW's can be present. We write as 

Ai+ = |Ai+|e*<\ (3.16) 

and introduce the phases Qj and c/yj defined by 

=ae, (3.17) 

According to the analysis of the free energy [Sec. |III C| 
and Appendix \uja\ = |A^+/ cos(0j)| is independent 
of a. We therefore obtain 

(5_(r)) = Yl 2|A^++| cos(Q, . r + e.Oe"'^^' 

3 

x[l + (-l)"tan(6lj)]. (3.18) 
^From Eq. ( 3.18[ ), we deduce 



(5,(r)) = Y 2|A^++| cos(</.j) cos(Q, . r + 9,) 

i 

x[l + (-l)«tan(0,)], 
{Sy{r)) = Y 2|A^++| sin(0,) cos(Q, . r + 6,) 

3 

x[l + (-l)"tan(0j)]. 



(3.19) 



The phase 4>j determines the polarization of the SDW's, 
while Qj gives their positions with respect to the under- 
lying crystal lattice. The free energy is independent of 
and Qj. 

The ratio of the SDW amplitudes on even and odd 
chains is given by the factor 5j = (1 + ta.n{9j))/{l — 
tan(6'j)) shown in Fig. ^. S remains close to one when V 
is weak {V < Vd). S 1 or 6 ^ 1 when V is strong 
(tan(6') ^ ±1 for oo) showing that the SDW's be- 
come mostly localized on even or odd chains (depending 
on the sign of 9). As expected, 0+ and 0_ have opposite 
signs so that S+ ~ 6Z^ ■ 

Consider first the SDW phase occuring for a weak an- 
ion potential V (interband pairing). There is a single 
SDW with a wave vector Q = {2kp, n/b). The spin mod- 
ulations on two neighboring chains are out-of-phase and 
the ratio 6 of their amplitudes is close to one (Fig. 10 1). 



The SDW's corresponding to intraband pairing are 
also commensurate in the transverse direction, but with 
{Qy)j = 7r/26: the spin modulations on two neighbor- 
ing chains are in phase quadrature. Fig. |lo| b shows the 
spin modulation in the phase with a single SDW at wave 
vector Q_. The spin modulations below T^°'^^ when two 
SDW's coexist is shown in Fig. p^ . The beating phe- 
nomenon is due to the fact that {Qx)- 7^ iQx)+- 



IV. CONCLUSION 



Anion ordering in the organic conductor 
(TMTSF)2C104 is expected to strongly influence the 
SDW instability. While a weak anion potential {V ^ tt) 
suppresses the SDW instability, a strong anion poten- 
tial leads to a rich phase diagram. When V becomes of 
the order of the interchain hopping amplitude (V ~ tb), 
the effective hopping between even and odd chains is 
reduced. This opens up the possibility to have two suc- 
cessive instabilities when the temperature decreases. The 
first one primarily occurs on even (or odd) chains and 
destroys the Fermi surface of one of the two electronic 
bands. The other electronic band remains metallic in 
this phase. At lower temperature, a second transition 
occurs, primarily on odd (or even) chains, making the 
whole Fermi surface gapped. 

According to Ref. |l^, the gap due to the anion poten- 
tial in the electronic dispersion is of the order of the trans- 
verse bandwidth, and the intraband nesting is almost 
perfect. In our model, this corresponds to a large value of 
V, i.e. V ^ tb. In this situation, one expects the ground- 
state to be a SDW. Experimentally, the ground-state is 
found to be a superconductor. There is no contradiction 
with our model though, since we did not consider the 
possibility of a superconducting transition. However, the 
observation of FISDW phases in (TMTSF)2C104 does 
require deviations from perfect nesting (otherwise, only 
the phase = would be observed) . This suggests that 
even when the anion potential is strong (large V in our 
model), important deviations form perfect nesting per- 
sist. The latter can be taken into account by considering 
tib in our model. Alternatively, the observation of the 
FISDW phases could indicate that the anion potential 
is weaker than what is predicted by quantum chemistry 
calculation. In our model, this would correspond to the 
region Vd < V < Vc2, where the ground-state is metallic. 
[Note that Vd = if t2b is large enough.] 

Two successive transitions have beep, observed in 
(TMTSF)2C104 at high magnetic field.Em The possibil- 
ity that each of these transitions corresponds to an in- 
traband pairing, thus partially destepying the Fermi sur- 
face, has been suggested early on.B Although our con- 
clusions do not apply stricto sensus to the experimental 
situation, since they are restricted to the zero-field case, 
they indicate that this scenario could indeed take place 
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in (TMTSF)2C104f4iugh field, in agreement with the 
conclusions of Refs-EHlBB 

Beside the existence of two successive transitions at 
low temperature when V is strong enough, an important 
result obtained in this paper is the overall phase diagram 
as a function of V. In particular, we have obtained a 
SDW-metal-SDW transition at low temperature (Figs. ^ 
and @). 

Finally, we note that our results seriously call into 
question the validity of most of the previous works 
on the SDW transition at high magnetic field in 
(TMTSF)2C104. In the presence of anion ordering, the 
SDW order parameter necessarily has two Fourier compo- 
nents, and the standard Stoner criterion cannot be used 
anymore. Instead, one should c onsid er th e gen eralized 
Stoner criterion obtained in Sec. [IB [Eq. ( ^.23 )]. 



Note added: After completion of this work, we bp=. 
came aware of a related work by Zanchi and Bjclis.E2l 
These authors have considered the effect of anion or- 
dering on the SDW instabihty in (TMTSF)2C104 and 
obtained the transition temperature within RPA. Their 
results are similar to ours. However, they did not take 
into account a finite t^i, and did not study the coexistence 
of the two SDW's with wave vectors Q-|- and Q_ at low 
temperature. 
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APPENDIX A: QUARTIC CONTRIBUTION F4 TO THE FREE ENERGY 



Using Eqs. (2^), ( 2.11 ) and (3J), the quartic part of the free energy [Eq. (3^)] can be written as in Eq. ( 3.1C ) with 



BaUl,j2,j3,j4;Pl,P2,P3,P4:) = f^^^ ^ ^ Fa{jl ■ ■ ■ js, Pl ■ ■ ■ Ps, P5 ■ ■ ■ p's', ky) ^aijl 



2N 



ky j5---j8 P5---PS p'^---p'f^ 

00 

^ 51/ ^a{pl-p2+P3+P■l)%- + (-P5-P6-P7-P8+p'^+p'(i+p'^+Ps)^■■n^' 
?l— — 00 

Aa(jl • • • JS; fcy) = ^ E - <k)"'(*^ - <k-aQ,,)"' 

X - <k-aQ,,+aQ,3 )"'(*'^ " «a,k-aQ,, 

'ky -aQj^ y+aQi3y-(a(P2 +P3 )-p'^ +P6 -p'g +P7) ^ 
' ky-aQj^y+aQj^y-{a{p2+P3)~p'^+P6-p'g +P7)-|- 

'ky-aQj-^y-{a(pi-l)-p5+p'g)§- 'fcj; -aQjii; - (q(pi - 1) -P5+Pg) ' 

where iV_L = Ly/b is the total number of chains. Since Vja — 0, we have 



A = c u- 

cos{9j), 



Noting that the conditio n Q, + Qj^ = Qj^ + Qj^ implies ji = j2 = js = ji or ji = ,72 = — js 
~ii = iii we obtain Eq. ( |3.11 ) with 



r,ky) 



-Ji or ji 



E "^JPi ^IP^ ^JP3 ^IP-i Ba{jJ,jJ;Pl,P2,P3,P4), 



Pl---Pi 



3 Pi ■■■Pi 

+ Cjp^ C-jp2 C-jp3 



Cjp4 Ba{j,-j,-j,j;Pl,P2,P3,P4)]- 



(Al) 



(A2) 



(A3) 

-.h = 



(A4) 



The main contribution to Bj comes from the type of diagrams shown in Fig. |^. This diagram involves the quantity 
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Kaiii = ■ ■ ■ = js = i;ky) ^ — ^ {iu 



a,k— aQj > 



(A5) 



To evaluate the order of magnitude of this diagram, we can assume perfect nesting. We then find 



Js 



1 

2^2 ■ 



A typical contribution to C is shown in Fig. ^. The order of magnitude of the diagram is given by 



Q,k— aQj 



a,k— aQj+aQ- 



(A6) 



(A7) 



To evaluate the preceding equation, we can assume perfect nesting in the j band: 

T 



Ka (j, j, - j, - j, j, j, h , j ; ky ) 



(A8) 



where |a ^ ^ ^ 
conclude that C/Bj = 0{T^/V^ 



V. We therefore 



APPENDIX B: MINIMUM OF THE FREE 
ENERGY F 

The minimum of the free energy is obtained by solving 
the equations dF/du*^ = 0: 



,^[A+ + 2B+\u+X + C\u_^\^]^0, 



2] =0. 



(Bl) 



In the metallic phase, Uja = and F = Q. Below 
T~ , there is a phase with it_Q, ^ and u^a ~ 0. From 
Eqs. @ and (|3l|), we deduce 



2BZ 



F_ = -- 



2K 



(B2) 



Let us now consider a phase wit h tw o coexisting order 
parameters: u^a,U-a 7^ 0. Eqs. (Bl) and (|3.13 ) yield 



\U+c 



-2A+B^ +A^C 
AB+B^ - ' 

-2A^B+ + A+C 
AB+B^ - C2 ' 



(B3) 



This solution is allowed only if \uja\ > 0. The phase with 
two coexisting order parameters will be observed only if 
it has a lower free energy than the phase with a single 
order parameter, i.e. i^coex 

-AiB- - AlB 



AB+B. 



C2 



A+A-C A^ 
— < — " 



2B. 



(B4) 



For AB^B^ — < 0, this condition can be rewritten as 
— (2A_Bj^ — A_|_C)2 > 0, which shows that coexistence 
is not possible. For AB^B^ — > 0, the condition 
-Fcoox < becomes -{2A^B+ - A+Cf < and is 
therefore always satisfied. We thus conclude that coexis- 
tence occurs when the condition 



AB+B^ - > 



(B5) 



is fulfilled. The corresponding transition temperature 
T^°'^'' is then determined from \uja\ > [Eqs. @]. 
When C/Bj < 1, T^^°'=^ is obtained from > 0. 

Writing Aj = aj{T — T-!) and neglecting the temperature 
dependence of Bj and C, we obtain 



rpCOeX rp-\ 



t: 



(B6) 



where 7 
0{T^/V^). 



a-C/{2a+B^) is a constant of order 
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FIG. 1. Possible nesting vectors in presence of the anion potential. Qintor (Q±) corresponds to interband (intraband) pairing. 
The magnitude of the band splitting at ky = ±7r/26 equals 21/, where V is the strength of the anion potential (see text). 
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FIG. 2. Susceptibility X°(c1iQ) in the absence of anion ordering (V = 0). Here x° = X+^ = X+ii a-nd qx is measured from 
2fcF-The maximum is located at the interband nesting vector Qinter with ((5inter)i/ ~ n/b. The susceptibilities shown in Figs. 
are obtained for T = 0.024;,, t2b = 0.1ft, and t^i, = 0.02f5. 





FIG. 4. Diagonal component x^CQi Q) of the susceptibility for V/tb = 1. Besides the maximum at Qinter, there are two peaks 
at Qy = 7r/2b corresponding to the intraband nesting vectors Q+ and Q-. The notations are the same as in Fig. 0. 




FIG. 5. Feynman diagrams for x within RPA. 
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1.5 



t2b= tsb = 



_ 1.0 




FIG. 6. Phase diagram with t2b = t^b = 0. The thick vertical line is a guide to the eye for estimating the critical potential at 
which the transition from Qinter to Q± takes pl ace. F or V > Vc — l.lib, there is coexistence of two SDW's (with wave vectors 
Q+ and Q-) below the transition line (see Sec. [II C). 
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Q inter 
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FIG. 7. Phase diagram with non-zero t2b and t^b showing the intermediate metallic region. Th e soh d line indicates the 
transition temperature when tsb ~ 0. Below y™™ 



Tc , the two intraband SDW's coexist (see Sec. [II C ) 




FIG. 8. Ratio S = {1 + tan(^))/(l — tan(S)) of the SDW amplitudes on even and odd chains (see Sec. [IID). The figure 
shows 5 if 5 < 1 and S^^ if 5 > 1. The parameters of the plot are the same as in Fig. ^ (with t-jb = 0.02t6). 
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FIG. 9. (a) Main contribution to Bj. (b) A typical diagram contributing to C. 




FIG. 10. Schematic representation of the spin modulation in the SDW phases (n denotes the chain index), (a) SDW 
phase for a weak anion potential V (interband pairing), (b) SDW phase for a strong anion potential (intraband pairing), (c) 
Low-temperature SDW phase (T < T™"") where two SDW's coexist. 
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